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Although canine respiratory coronavirus (CRCoV) is an important respiratory pathogen that is prevalent in many 
countries, only one complete genome sequence of CRCoV (South Korea strain K37) has been obtained to date. 
Genome-wide analyses and recombination have rarely been conducted, as small numbers of samples and limited 
genomic characterization have previously prevented further analyses. Herein, we report a unique CRCoV strain, 
denoted strain BJ232, derived from a CRCoV-positive dog with a mild respiratory infection. Phylogenetic 
analysis based on complete genome of all available coronaviruses consistently show that CRCoV BJ232 is most 
closely related to human coronavirus OC43 (HCoV-OC43) and BCoV, forming a separate clade that split off early 
from other Betacoronavirus 1. Based on the phylogenetic and SimPlot analysis we propose that CRCoV-K37 was 
derived from genetic recombination between CRCoV-BJ232 and BCoV. In detail, spike (S) gene of CRCoV-K37 
clustered with CRCoV-BJ232. However orflab, membrane (M) and nucleocapsid (N) genes were more related to 
Bovine coronavirus (BCoV) than CRCoV-B232. Molecular epidemic analysis confirmed the prevalence of CRCoV- 
BJ232 lineage around the world for a long time. Recombinant events among Betacoronavirus 1 may have 
implications for CRCoV transmissibility. All these findings provide further information regarding the origin of 
CRCoV. 


1. Introduction 

Canine respiratory coronavirus (CRCoV), an enveloped positive- 
stranded RNA virus, belongs to the betacoronavirus genus, subgroup 2a 
(Erles et al., 2007; Erles et al., 2003). Unlike the enteric canine 
coronavirus (CCoV) belonging to the alphacoronavirus genus, CRCoV 
is associated with mild to severe respiratory signs and proposed as an 
etiological agent of canine infectious respiratory disease (CIRD) 
(Buonavoglia and Martella, 2007). CRCoV was first detected in 2003 
by RT-PCR in the UK from dogs suffering from severe respiratory 
disease (Erles et al., 2003). Later, the CRCoV-4182 strain was success¬ 
fully isolated from a dog with respiratory signs and was cultured in the 
HRT-18 cell line (Erles et al., 2007). To date, data have been made 
publicly available for only one complete genomic sequence (CRCoV- 
K37 strain), which was isolated from a South Korea dog (An et al., 
2010; Lim et al., 2013). 


CRCoV was included in species Betacoronavirus 1, which also 
contain BCoV, equine coronavirus, HCoV-OC43, porcine hemaggluti- 
nating encephalomyelitisvirus (PHEV), human enteric coronavirus 
(HECoV) and CRCoV (King, 2011). In contract, Canine coronavirus 
(CCoV) is a member of alphacoronavirus, which include feline corona¬ 
virus and porcine transmissible gastroenteritis among others. The virus 
was distinct from CCoV and only showed 69% nucleotide identity in 
polymerase and 21% amino acid identity in spike protein (Buonavoglia 
and Martella, 2007; Erles and Brownlie, 2008). 

The evolutionary origin of CRCoV remains uncertain, although 
successful experimental infection of puppies with BCoV has suggested 
that CRCoV is of bovine origin and transmitted to dogs from cattle 
(Kaneshima et al., 2007 ). Meanwhile, HCoV-OC43 may have emerged 
after viral transmission from cattle to people (/ijgen et al., 2005). The 
high nucleotide identity between CRCoV and BCoV, HCoV-OC43 and 
HECoV indicated that the four viruses have a common ancestor ( irles 
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and Brownlie, 2008), and demonstrated the occurrence of repeated 
host-species shifts ( Tjgen et al., 2006). Although CRCoV strains are 
frequently circulating in farmed and pet dogs worldwide (Erles and 

Brownlie, 2005; Kaneshima et al., 2006; Priestnall et al., 2006; 
Priestnall et al., 2007; Yachi and Mochizuki, 2006), causing respiratory 
disease, little is known about their genomic evolution and recombina¬ 
tion. In the present study, we sequenced and analyzed the first complete 
genome of CRCoV derived from China (CRCoV-BJ232 strain) and 
discuss its role in CRCoV evolution. 

2. Materials and methods 

2.1. Ethics statement 

This research was reviewed and approved by the ethics committee 
of the National Institute for Viral Disease Control and Prevention of the 
Chinese CDC. All animals were treated strictly according to the guide¬ 
lines for the Laboratory Animal Use and Care from the Chinese CDC. 

2.2. Specimen collection and handling 

From December 2013-March 2014, 246 swab specimens (nasophar¬ 
yngeal and anal swabs) were collected from diseased dogs at the Animal 
Hospital of China Agricultural University. Most of the dogs showed 
symptoms of gastrointestinal and/or respiratory tract infections, in¬ 
cluding coughing, lacrimation, vomiting and diarrhea. All the swab 
samples were collected by professional veterinarians and were stored at 
-80 °C. 

2.3. CRCoV detection 

For molecular detection, total RNA was extracted from swab speci¬ 
mens using a QIAamp Viral RNA Mini Kit (QIAGEN, Germany) 
according to the manufacturer’s instructions. The RNA was used as a 
template for reverse transcription-PCR (RT-PCR), and the partial 
sequences of the S, RNA-dependent RNA polymerase (RdRp), and N 
genes were amplified using gene-specific primers (Table SI). In detail, 
reverse transcription was performed by using a Superscript III kit 
(Invitrogen, San Diego, CA). The PCR mixture (25 pi) contain 1 pi 
cDNA, 1 pi forward primer(10 pM), 1 pi reverse primer (10 pM), 12.5 pi 
Q5 @ High-Fidelity 2X Master Mix (New England Biolabs) and 9.5 pi 
nuclease-free water; the mixtures were amplified with 30 cycles at 98 °C 
for 10s, 52 °C for 30s, and 72 °C for 1 min and final extension step at 
72 °C for 2 min in an automated thermal cycler (BIO-RAD, T100 
Thermal Cycler). 

2.4. Complete genome sequencing of CRCoV 

To better characterize the CRCoV discovered here, the CRCoV- 
positive samples were pretreated and used for random amplification as 
described previously (Wang et al., 2016); the amplified DNA was used 
as a template for Illumina HiSeq2500 sequencing (2 x 125 bp, paired- 
end reads). The sequencing procedure was conducted by Beijing Berry 
Genomic Company. The full-length genome of CRCoV was assembled 
using SOAP de novo and was annotated based on the sole annotation of 
CRCoV K37 strain (JX860640.1). 

2.5. Phylogenetic analysis 

To determine the evolutionary relationships between the CRCoV 
discovered here and those identified previously, we collected most of 
the complete sequences of the Betacoronavirus 1 and conducted 
nucleotide acid sequence alignments using the MAFFT algorithm. The 
phylogenetic analyses were conducted based on the nucleotide acid 
sequences of the orflab, S, M, and N genes and complete genome of 
CRCoV. All the phylogenetic trees were constructed in the PhyML 


software using the maximum likelihood method, with 500 bootstrap 
replicates. The following Betacoronavirus 1 genome sequences from the 
GenBank database were used for the phylogenetic analysis as follows: 
BCoV reference strain (NC_003045), HECoV (NC_012950), HCoV-OC43 
(NC_005147), CRCoV-37 (JX860640), and BCoV strain Quebec 
(AF220295). 

2.6. Recombinant analysis 

To confirm the potential recombination events, similarity plots and 
bootscanning analyses were conducted using the SimPlot software; a 
sliding window of 1000 nucleotides and 100-nucleotide steps were used 
as the default settings. In addition, a schematic diagram of the complete 
CRCoV genome was examined to determine the recombination sites. 

2.7. Nucleotide sequence accession number 

The complete genome of ChinaBJ232 was deposited in the GenBank 
database under accession no. KX432213. Other sequences generated in 
this study, including the partial S, RdRp and complete N gene 
sequences, were also assigned accession numbers: KX432214- 
KX432217. 

3. Results 

3.1. Identification of CRCoV 

A total of 246 swab samples were collected in Beijing from 
December 2013 to March 2014.RT-PCR was performed to detect 
CRCoV, and 16 (6.5%) out of 246 swabs originating from Chinese dogs 
tested positive for CRCoV. In addition, five CRCoV-positive samples 
with high titers were initially processed for complete genome sequen¬ 
cing using deep sequencing technology. 

3.2. Complete genomic characterization of the BJ232 strain 

Due to the limited specimen size and low viral load, five clinical 
samples were used for deep sequencing, and only one complete 
genomic sequence was eventually obtained. The genomic data, formally 
referred to as “Canine Respiratory Coronavirus strain BJ232” and in 
this report as “CRCoV-BJ232,” were deposited in GenBank (KX432213). 
The CRCoV-BJ232 genome isolated in China was found to have 30,868 
nucleotides and showed a genomic organization that was highly similar 
to the Betacoronavirus 1 species: 5'-UTR(nt l-143),ORFlab (nt 
144-21427), ORF2 (nt 21437-22273), HE gene (nt 22285-23559), S 
gene (nt 23574-27665),E gene (nt 28355-28609), M gene (nt 
28624-29316), N gene (nt 29326-30672), and 3'-UTR (nt 
30673-30868). Between the S and E gene, three small ORFs were 
identified: ORF4a (nt 27655-27789), ORF4b (nt 27825-27902) and 
ORF4c (nt 28039-28293), encoding 4.9-, 2.7- and 12.8-kDa non- 
structural proteins, respectively, which were analogous to those noted 
in strain CRCoV-K37. Meanwhile, conventional PCR and Sanger 
sequencing were conducted to verify the complete genome sequence 
of CRCoV-BJ232 with 100% nucleotide identity. Thus far, there are 
only two CRCoV strains with full-length genomic information submitted 
to GenBank; the previously reported CRCoV-K37 and our current strain, 
CRCoV-BJ232. Unfortunately, attempts to isolate CRCoV-BJ232 from 
the five positive samples were unsuccessful; future studies to isolate the 
virus from more canine samples will allow characterization of its 
pathogenicity and transmissibility. 

3.3. Phylogenetic analysis 

Based on the analysis of whole-genome phylogeny, CRCoV-BJ232 
was classified into Betacoronavirus 1 species of betacoronavirus group 
2a ( 7 ig. 1A) and shares a common origin with BCoV, CRCoV-K37, 
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Fig. 1. Phylogenetic analysis of all available coronaviruses with reference sequence based on complete genome. A, the current taxonomy of the family Coronaviridae; B, the evolution of 
Betacoronavirus 1 based on complete genomes. StrainBJ232 (accession no. KX432213, indicated by the red dot) was the first CRCoV with complete genome derived from China. The 
previously reported CRCoV-K37 strain, which was detected in South Korea, was the sole other CRCoV strain with complete genome sequence available. Some other representative strains 
of the Betacoronavirus genus, subgroup 2a, were also analyzed in PhyML program using the maximum likelihood method, with 500 bootstrap replicates. The numbers at the nodes 
represent the bootstrap support. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) 


HECoV and HCoV-OC43. Moreover, the latter four viruses clustered 
together with high branch bootstrap values ( Fig. IB). To investigate the 
genetic relationships among the Betacoronavirus 1 species, we gener¬ 
ated phylogenetic trees based on the sequences oforflab, S, M and N 
genes ( ug. 2). Phylogenetic analysis based on different genes consis¬ 
tently show that CRCoV BJ232 is most closely related to HCoV-OC43, 
forming a separate clade of CRCoV that split off early from other 
Betacoronavirus 1. However, CRCoV-K37 strain closely clustered with 
CRCoV-BJ232 strain in S gene, and was more related to BCoV strain 
Quebec (AF220295.1) in the region of orflab, M and N genes. The 
major genes of the CRCoV-BJ232 strain showed a consistent trend in 
the evolution, while the phylogenetic topology conflicted in the sub- 
regions of CRCoV strain K37, especially in the S and orflab genes. 
These findings are consistent with recombination, a phenomenon not 
uncommon in coronaviruses, implying that the CRCoV K37 strain may 
have arisen from genetic recombination. 

3.4. Recombinant analysis 

Similarity plots and bootscanning analyses ( ’ig. 3) were performed 
to confirm the recombination between CRCoV and BCoV. The pre¬ 
viously reported CRCoV K37 strain was used as a query sequence, and 
was compared with HCoV-OC43 (NC.005147), BCoV (BCoV- 
RefNC_003045, BCoV-DB2 DQ811784, BCoV-Quebec AF220295), HE¬ 
CoV (HECoV-Ref NC_012950, HECoV-4408 FJ415324), and CRCoV 
(CRCoV-BJ232 KX432213). In the S gene, strain K37 possessed the 
highest similarity with the BJ232 strain between position 23,842 and 
29,072, as expected. However, in the orflab, HE, M and N genes, strain 


K37 was apparently related to BCoV, which further confirmed the 
occurrence of genetic recombination between CRCoV and BCoV. The 
consistency in the results of bootscanning and phylogenetic analyses 
also supported the possibility of recombination in the S gene. 

3.5. Prevalence of recombinant CRCoV worldwide 

We collected all available CRCoV sequences from the GenBank 
database and conducted comprehensive phylogenetic analyses sepa¬ 
rately based on HE, S, M and N genes to understand the epidemiology 
distribution of the recombinant strain(South Korea strain K37). The 
BJ232, BJ221, BJ34 and BJ202 strains were collected from China in 
this study, while the K37, K39 and K9 strains came from South Korea, 
and the 240/05 and 4182 strains were derived from Italy and the 
United Kingdom, respectively. As shown in Fig. 4, the CRCoV strains 
from China, including the BJ232, BJ221, BJ34 and BJ202 strains, 
clustered with strains prevalent in Italy and the United Kingdom (the 
240/05 and 4182 strains), and the CRCoV strains detected in South 
Korea (K37, K39 and K9 strains) were all derived from recombinant 
lineage. 

4. Discussion 

CRCoV is a causative agent for CIRD and a novel pathogen detected 
in respiratory samples from dogs (Erles and Brownlie, 2008). A 
serological survey of CRCoV conducted worldwide revealed that 
antibodies against CRCoV were common in the canine population 

(Erles et al., 2004; Isaacs et al., 1983; Pratelli et al., 2003; Priestnall 
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Fig. 2. Evolutionary analyses of major coding genes of Betacoronavirus 1. The coding genes used in the phylogenetic analyses included orflab, S, M and N. The two CRCoV strains (BJ232 
and K37 strains) are indicated in red. All the sequences were analyzed in PhyML program using the maximum likelihood method, with 500 bootstrap replicates. 


et al., 2006). In contrast, corresponding genome-wide analyses have 
rarely been conducted (An et al., 2010; Lorusso et al., 2009). In the 

present study, we characterized the first complete genome of CRCoV 
(named strain BJ232) derived from China, and phylogenetic analyses 
showed that CRCoV-BJ232occupied a deep branch at the root of 
members of Betacoronavirus 1, being distinct from HCoV-OC43 and 
BCoV. Another previously described CRCoV strain, K37, was identified 
to have arisen from recombination between CRCoV-BJ232 and BCoV. In 
detail, the CRCoV-K37 strain was related to CRCoV-BJ232 in the S gene, 
while in the orflab, HE, M, and N genes, CRCoV-K37 was noted to 
possess a close relationship with BCoV. 

Recombination is a common phenomenon in coronaviruses and 
thought to contribute to the emergence of new pathotypes (Gorbalenya, 
2008; Wang et al., 2015). Thus far, most of the recombination events of 
coronaviruses have been reported between species of the same group 
(Herrewegh et al., 1998; Keck et al., 1988), such as among the bat- 
associated CoVs (Corman et al., 2014) and 229E-related CoVs (Corman 
et al., 2015), and the major recombination breakpoint has mainly been 
within the S gene. Genetic recombinant among alphacoronavirus 1 
species frequently occurred, such as ferret coronaviruses (Lamers et al., 


2016; Minami et al., 2016). However, few researches about genetic 
recombination of Betacoronavirus 1 species have been conducted. The 
discovery of genetic recombination between CRCoV and BCoV confirms 
the phenomenon, and intraspecies recombination therefore exists for 
Betacoronavirus 1 species. 

Phylogenetic analyses based on all available CRCoV genes showed 
that the two clusters of CRCoV, including recombinant strain (South 
Korea K37 strain) and the Chinese BJ232 strain, have been prevalent 
worldwide for many years. CRCoV and BCoV were closely related in 
Betacoronavirus 1, located on the outer leaves of phylogenetic tree, and 
form an independent branch. It was difficult to determine the evolution 
orientation (transmission from cattle to canine or both from murine), 
although the latest research supports the murine origins of 
Betacoronavirus 1 (Lau et al., 2011, 2015). More CRCoV genome 
sequences are needed to understand the evolutionary relationship 
between CRCoV and BCoV. 

Regarding the pathogenicity and transmissibility of CRCoV, the 
Chinese BJ232 strain and South Korea K37 strain were both derived 
from diseased dogs with respiratory infections. Genomic analyses 
showed that the S genes, which mediate coronavirus attachment, were 
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Fig. 4. Phylogenetic analyses of all available CRCoV strains based on HE, M, N and S genes. The BJ232, BJ221, BJ34 and BJ202 strains come from China in this research,K37, K39 and K9 
strains were collected from South Korea, and the 240/05 and 4182 strains were derived from Italy and the United Kingdom, respectively. All the genes were analyzed by PhyML program 
using the maximum likelihood method, with 500 bootstrap replicates. 


highly similar in these strains. Because recombination may be a driving 
force for the formation of pathogenic viruses from less pathogenic virus 
("orni et al., 201' ), recombinant may play an important role in the 
pathogenicity and transmissibility of CRCoV. Limitation of this study 
was that the isolation of CRCoV-BJ232 failed on cell culture. More 
CRCoV strains were needed to be isolated or rescued using reverse 
genetic technique to perform experimental inoculation of dogs. 

In conclusion, herein, we reported the first complete genome 
sequence of a CRCoV derived from a diseased dog in China, and 
confirm that intraspecies recombinant exists for Betacoronavirus 1 
species. The results may have implications for the origin and evolution 
of CRCoV. 
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